Vector potential analysis of the helicon antenna in vacuum 
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The helicon antenna is a well-known device in the field of electric propulsion. Here we investigate 
the vector potential produced in vacuum by such an antenna with typical size parameters. Both 
a static and a dynamic analysis are performed. The dynamic calculation is evaluated at both the 
usual operating frequency and one which is slightly greater. At the higher frequency, a pulse of 
electromagnetic energy is found to propagate along the cylindrical axis in either direction. The 
possible adaptation of the helicon antenna as a RF injection device for burning plasma is discussed. 

PACS numbers: 52.25.Jm, 52.50.Dg, 52.40.Fd 



I. INTRODUCTION 

The helicon antenna, consisting of two rings separated 
and joined by two curving limbs, is a well-known device in 
its application as an ion source for electric propulsion fl^ 
Driven by a radio- frequency amplifier, its energy cou- 
ples to the plasma medium to drive the dissociation of 
propellant material, such as argon gas. Here we inves- 
tigate the vector potential produced by such a device in 
vacuum, with an eye towards its adaptation as a source 
of RF energy for injection into a burning plasma. 

Our methodology is based on the numerical evalua- 
tion of solutions to the inhomogeneous field equations 
— —fioJ'^ in Lorenz gauge d^A^^ = 0, where 
= (a/co,A) is the four- vector for the potential, 
Jf^ = (cqJ, J) is the four-vector for the current source, 
and Cq^ = yuoeo gives the speed of light in vacuum. The 



d'Alembertian operator is □ = d'^di, 



+ in 



a frame with metric signature ( — ,-|-, -|-). We subscribe 
to the reductionist viewpoint [a, 01 that the potential is a 
more fundamental description than the field formulation, 
and that the Lorenz gauge makes explicit the relation be- 
tween the continuity of the potential and the continuity 
of the source 9^ = 0. 

The solution for a given J** is found by directly 
imaging the current source throughout the region, by 
which we mean that the contribution from each source 
element J^{t' ,v') is calculated at a point (t, r) with re- 
gard to the propagation delay, or retarded time t' — 
t — Ar/ca, as well as the spatial distance to the element 
Ar ~ |r — r'|. The antenna is assumed to be driven by a 
feed current arranged so that its contribution is negligi- 
ble, leaving to calculate only the potential produced by 
the antenna. We neglect any charge accumulation along 
the conductor, so that the static charge density vanishes, 
thus J = and V • J — 0. The solution to the wave 
equation for a line source is then written 



Ac(i,r) 



47r 



Ic(i',r') 



dl 



(1) 
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where the subscript C reminds us that these vec- 
tors must be expressed in Cartesian components A — 
{Ax,Ay,Az), neglecting any self-interaction between 
different parts of the antenna. 

The antenna is modeled as a collection of discrete line 
current elements Ifc = I(r'j,), so that the integral becomes 
a sum J^i dl J2k ^k, where Afc is the length of the fcth 
element. The antenna is assigned typical size parameters 
of half- width w = OA m and half- height /i = 0.2 m (or 
1 and 2 dm respectively), and the feed current is nor- 
malized to iiqIq/Att = 1. We assume that the driving 
current is provided along a twisted pair cable so that the 
opposing phases produce canceling contributions to the 
overall potential. The Cartesian {X, Y, Z) and cylindri- 
cal {R, 4>, Z) coordinates are related by X = i? cos and 
Y = Rs'mcj). The static analysis calculates the poten- 
tial in the XY plane at various Z along the cylindrical 
axis. The dynamic analysis calculates the potential in 
the ZR plane at various angles around the cylindrical 
axis. The propagation analysis computes A in the ZR 
plane at cj) — -k /A and various times t. Our methodology 
is a straightforward implementation of classical electro- 
dynamics in the potential formulation, such as found in 
Ref. B. 



II. STATIC ANALYSIS 

For the static, time independent analysis, the Lorenz 
gauge is equivalent to the Coulomb gauge, V ■ A = 
0. Position r'(A, 0') along the helicon antenna may be 
parametrized by the azimuth (p' and the line distance 
along the conductor A to the location where the driving 
current is attached, here taken to be at {X' ,Y\ Z') = 
(-u;,0, -h). The ring at Z' = h will be called the "head", 
and the one &% Z' = —h will be the "tail". The vertical 
limbs we call "arms" , and the geometry of the antenna is 
depicted in Fig. [1] where the "legs" carrying the driving 
current are not shown. It takes two revolutions around 
the circle to cover the entire antenna, thus the positive 
half is identified with 0'^ S [0, 27r] and the negative with 
€ [0, —2-k]. The azimuth is parametrized over iV^ val- 
ues, so that = 2-K/N(f,. The line elements for the rings 
have a length = wA^, and those for the arms have a 
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v = MHz 



III. DYNAMIC ANALYSIS 
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FIG. 1. The helicon antenna is modeled as a collection of 
discrete line current elements as described in the text. The 
geometry is shown in units of 1 dm = 0.1 m throughout. 



length Afe = (A 



w 



U/2 



where Az — 2hA^/TT is 



the vertical extent for arm elements at Z'f. — 2h<p'^./Tr — h. 
Alternately, one could take J2k ^ E^' Ehcad,taii,arm 
one revolution of the circle (j)' G [0, Stt]. 

By Kirchhoff 's point rule, the magnitude of the cur- 
rent along either branch of each ring is Ir = /o/2, and 
along the arms it is la = lo- The arms must cover a 
height of 2h over an azimuthal range of tt, thus the ratios 
Iz/I<p = 2h/TTw and {Ia/14,)'^ = 1 + {Iz/I<p)^ determine 



I^ 



-Id, sin (h'l, and /y 



cos for /, 



(2ft/7rw)^]~^/^ so that la /o(- sin0'^, cos 0'^, 2/1/™), 
and similarly for I,, with sign changes as appropriate for 
the two branches of each ring. The contribution from 
each element I(rjj) to the solution A(r) is calculated by 
multiplying the current 1^ by Ak / Ar for each component 
in {X,Y,Z). 

On the midplane of the antenna Z = dm, the vector 
potential A has its greatest magnitude at the locations 
of the arm elements (0, ±w, 0), as seen in Fig. [2j and is 
dominated by the X and Z components as expected; the 
Y component is about a tenth of the magnitude of the 
others. At three quarters of the antenna height, Z = 1 
dm shown in Fig. |3l the locations of the arms are clearly 
discerned in the vector potential. In the plane of the head 
ring, Z = 2 dm shown in Fig. |31 the vector potential is 
mostly in the azimuthal direction, with a strong vertical 
component only at the junctions of the head ring with 
the arms. Away from the antenna, at Z = 3 dm shown 
in Fig. [SJ a significant vector potential still exists, with 
both vertical and azimuthal components. 



Let us now consider driving the antenna with an os- 
cillatory current of frequency = l/r > 0, whose wave- 
length in vacuum is Aq = cq/v. The phase of the cur- 
rent along the antenna is determined by the line distance 
along the conductor A ((/)') from the current element to the 
location selected as the zero phase reference, where the 
feed current attaches to the antenna. The antenna has a 
half-length of A^ax = A(±27r) = ttw+Httw)'^ + {2/i)2]i/2. 
The way we have orientated our antenna, the posi- 
tive branch consists of the tail ring and the arm with 
1^ < 0, while the negative branch is the head ring and 
the arm with y > 0. To reduce the number of brack- 
ets, let us denote expi(Q;) = exp(za). The complex 
current along the positive branch of the conductor is 
I+{t,\) = /q expi[27r(A/Ao — i/r)] /o(t) expi(6'A), where 
7o(<) = loexpi{0t) for 9t = -27rt/T and Ox = 2ttX/Xo. 
Along the negative branch of the conductor, the cur- 
rent is /_(t, A) — —I+{t,\). Propagation delay in- 
troduces the retarded phase 9r — 27rAr/Ao such that 
I+{t',X) = Ioexpi{6x + Or)- Here, /q equals Ir or la 
as appropriate. The complex vector potential Afe(i,r) 
is calculated by imaging each discrete element I(i,r'j,) 
throughout the region with weighting Ak/Ar and net 
phase expi(6'fc + 9t + Or) for 0^ = 27rAfc/Ao, so that the 
physical potential is A = Re A^. 



A. Low frequency 

Considering a driving frequency oi v — 13.54 MHz 
commonly used to excite argon gas, the vacuum wave- 
length is Ao 22 m, and Amax/Ao = 3.7%. The first 
step is to calculate the phase of the current along the 
antenna Ok, which we show in Fig. [6] The current is 
then expressed in Cartesian coordinates I = {Ix,IytIz) 
for each element at rj. with phase Ok- Next one needs 
the distances Ar(r) throughout the region from r'^, to r. 
From Ar one gets the magnitude weighting A^/Ar and 
the retarded phase Or- We set t/r = so that we can 
investigate the spatial distribution of the vector potential 
in the ZR plane at various azimuths 0. 

Along the X axis, where </) = 0, the vector potential 
has the structure for Z,R > shown in Fig. \7\ The 
magnitude |A| is dominated by contributions in the Y 
and Z directions at the location of the junction of the arm 
and the head ring, with a component in the X direction 
for the interior of the antenna R < w- In the plane at 
(j) — tt/A in Fig. [51 the strongest contribution is along Y 
at the height of the head ring. On the Y axis, (j) = tt/2 in 
Fig. [HI the potential is mostly in the X and Z directions, 
with a component in the Y direction for the antenna 
interior. At (p = 37r/4 in Fig. (TUl the potential has a 
strong Z component from the antenna arm, with some 
additional structure along X and Y. 




Z= 1 Z = 3 

FIG. 3. Static vector potential at Z = 1. FIG. 5. Static vector potential sd Z — 3. 



v= 13.54 MHz 
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FIG. 6. For a driving frequency of u — 13.54 MHz, the phase 
difference at the maximum distance Amax from the location of 
the current feed legs is 3.7%. The colorbar indicates the phase 
in units of degrees, and the length of the arrow represents the 
magnitude of the real current at t/r = 0. 



B. High frequency 

Now let us suppose that the antenna is driven at a 
higher frequency, say v = 541.6 MHz, whose vacuum 
wavelength is considerably shorter Aq ^ 0.55 m so that 
Amai/Ao = 149%. The phase structure of the current 
along the antenna is displayed in Fig. 1111 where we see 
that the magnitude of the real current is zero where 9k + 
6t = ±7r/2. Repeating the azimuthal scan at the same 
angles as above for t/r = 0, the results are displayed 
in Fig. [1^ through Fig. [131 The vector potential has 
a complicated structure, the most interesting feature of 
which is the strong magnitude of A which develops at 
azimuths = 7r/4 and (j) — 37r/4. 

The complicated phase structure along the antenna in- 
dicates that we are in a regime in which interference ef- 
fects can lead to the propagation of an electromagnetic 
pulse. Restricting attention to the ZR plane at azimuth 
(j) = 7r/4, let us now investigate how the pulse evolves 
through time. Considering the set of times given by 
t/r G [l,8]/8, we display the results in Figs. [T6l through 
1231 respectively. Times separated by one half period have 
the same magnitude for A and vector components of op- 
posite sign. At the end of one cycle, the potential is back 
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(|)/7i=0/4,t/T=0/8 



FIG. 7. Vector potential for u = 13.54 MHz at = and 
t/r = 0. 



■., (a) 






(|,/7i= 1 /4,t/T=0/8 

FIG. 8. Vector potential for u = 13.54 MHz at cj) = 7r/4 and 
t/r = 0. 
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(|)/7i=2/4,t/T=0/8 



FIG. 9. Vector potential for v = 13.54 MHz at = 7r/2 and 
tjr = 0. 





(|,/7i=3/4,t/T=0/8 

FIG. 10. Vector potential for v = 13.54 MHz at = 37r/4 
and tjr = 0. 



V = 541.6 MHz 
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FIG. 11. For a driving frequency oiv = 541.6 MHz, the phase 
difference at the maximum distance Amax from the location of 
the current feed legs is 149%. The colorbar indicates the phase 
in units of degrees, and the length of the arrow represents the 
magnitude of the real current at tjr = 0. 



to the configuration found for 6'f = 0, in that Fig. [23] 
is the same as Fig. IS] What we find most interesting 
is how a pulse develops in the central region of the an- 
tenna, best seen in |A|, which then propagates along Z. 
The two peaks along the Z axis at i? = in Fig. [16] panel 
(d) represent the rectification of the peak and trough of 
the propagating pulse. The pulse is seen, by following its 
development through the figures, to detach from the an- 
tenna and continue into space beyond the head ring. An 
animation with twice the temporal resolution is available 
as an online supplement. 



IV. DISCUSSION 

The original motivation for this study was simply to 
develop an interesting exercise in classical electrodynamic 
field theory, with possible utility to investigators in elec- 
tric propulsion. However, by playing around with the 
driving frequency, we have found evidence that a signifi- 
cant pulse may be generated which propagates away from 
the antenna. We expect interesting results when the ratio 
Amax/ Ao is at half-integer multiples, so that constructive 
interference from the head and tail rings leads to a prop- 
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FIG. 12. Vector potential for u = 541.6 MHz At cj> = and FIG. 14. Vector potential for u = 541.6 MHz at </> = 7r/2 and 



t/r = 0. 



t/r = 0. 
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(|,/7i= 1 /4,t/T=0/8 



(|,/7i=3/4,t/T=0/8 



FIG. 13. Vector potential for u = 541.6 MHz at = 7r/4 and FIG. 15. Vector potential for ly = 541.6 MHz at = 3n/4 



t/r = 0. 



and t/r = 0. 
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FIG. 16. Vector potential for u = 541.6 MHz at = 7r/4 and FIG. 18. Vector potential for u = 541.6 MHz at (j> = n/4 and 



t/r = 1/8. 



t/r = 3/8. 
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(|,/7i= 1 /4,t/T = 2/8 



(|,/7i= 1 /4,t/T = 4/8 



FIG. 17. Vector potential for u = 541.6 MHz at = 7r/4 and FIG. 19. Vector potential for u = 541.6 MHz at </> = n/4 and 



t/r = 1/4. 



t/r = 1/2. 
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(|,/7i= 1 /4,t/T=5/8 



(|,/7i= 1 /4,t/T=7/8 



FIG. 20. Vector potential for v = 541.6 MHz at = 7r/4 and FIG. 22. Vector potential for = 541.6 MHz at t/i = 7r/4 and 



tjr = 5/8. 



t/r = 7/8. 







(|,/7i= 1 /4,t/T = 6/8 



(|,/7i= 1 /4,t/T=8/8 



FIG. 21. Vector potential for u = 541.6 MHz at = 7r/4 and 
t/T = 3/4. 



FIG. 23. Vector potential for u = 541.6 MHz &t <f> = tt/4 and 
t/r = 1. Online this figure is linked to a short video showing 
the temporal evolution of the vector potential. 
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agating pulse. Noticing the similarity between the higher 
frequency here and the range spanned by existing RF in- 
jection plasma heating technologies based on resonance 
with the ion cyclotron and the lower hybrid frequencies, 
together spanning a range of about 0.05 to 10 GHz 9], 
we wonder if the helicon antenna might be adapted for 
use in fusion engineering. For that purpose, the inclusion 
of a ground plane to achieve unidirectional propagation 
might be necessary to allow for mounting on the outside 
of the vacuum vessel, or perhaps the antenna unaltered 
can be placed towards the interior of the central solenoid 
in a tokamak so that both propagation pulses can be uti- 
lized. 

The model here has made several simplifying approxi- 
mations which need to be more fully developed. As the 
material of the conductor is unspecified, we have used 
the vacuum values of Eq and /iq in the determination of 
the phase along the antenna 0\. For greater veracity, 
the values of e and /i appropriate for the chosen material 
should be used instead. A volume or surface current den- 
sity could be used for the antenna rather than the line 
current used here, with straightforward modifications to 
the evaluation of the integrals. More interesting but more 
difficult is the incorporation of self interactions between 
different parts of the antenna — the current Ife cannot be 



immune to the effects of the current at other locations I^/ , 
which are beyond the scope of this article. Finally, one is 
interested in how the propagating electromagnetic pulse 
would interact with a medium of ionized gas located in 
the region beyond the antenna rings. By properly tuning 
the antenna dimensions and driving frequency, one hopes 
to drive resonances inside the plasma medium. 

V. CONCLUSION 

In this article we have presented a numerical evaluation 
of the vector potential produced by the helicon antenna 
in vacuum for both the static case and the dynamic case 
of nonzero driving frequency. We have considered both 
a typical driving frequency and one somewhat greater, 
where the higher frequency produces a significant elec- 
tromagnetic pulse which propagates along the cylindrical 
axis. The frequency of the propagating pulse for typical 
size parameters lies within the range spanned by ion cy- 
clotron and lower hybrid resonances, suggesting that an 
adaptation of the helicon antenna design might be use- 
ful for RF injection plasma heating. Extensions of the 
model for a more realistic description of an experimental 
apparatus are outlined. 
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